Counting statistics of interfering Bose-Einstein condensates 
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A method is presented that is able to predict the probability of outcomes of snapshot measure- 
ments, such as the images of the instantaneous particle density distribution in a quantum many-body 
system. It is shown that a gauge-like transformation of the phase of the many-body wave function 
allows one to construct a probability generating functional, the Fourier transform of which with 
respect to the "gauge" field returns the joint probability distribution to detect any given number 
of particles at various locations. The method is applied to the problem of interference of two inde- 
pendent clouds of Bose-Einstein condensates, where the initially separated clouds with fixed boson 
numbers expand and the density profile image of the overlapping clouds is registered. In the limit of 
large particle numbers, the probability to observe a particular image of the density profile is shown 
to be given by a sum of partial probability distributions, each of which corresponds to a noisy 
image of interference of two matter waves with definite phase difference. In agreement with earlier 
theoretical arguments, interference fringes are, therefore, expected in any single shot measurement, 
the fringe pattern randomly varying from run to run. These results conform to the physical picture 
where the Bose-Einstein clouds are in spontaneously symmetry broken states, the hidden phases of 
which are revealed by the density profile measurement via the position of the interference fringes. 



Shortly after the achievement of creating a Bose- 
Einstein condensate (BEC) in a dilute gas of ultra-cold 
atoms, an experiment was performed where two inde- 
pendent clouds of BEC were brought to overlap in space, 
and the particle distribution of atoms measured. A snap- 
shot picture of the density of the two overlapping clouds 
showed interference fringes analogous to the character- 
istic interference pattern of a double-slit experiment [l[ . 
The text-book double-slit case is a repeated measurement 
performed on the same single system, and the fringes in 
the arrival pattern of the matter wave emerges only af- 
ter many runs. On the contrary, the former interference 
fringes are present in a single snapshot picture of the 
density distribution of the overlapping condensates, with 
a phase, as determined from the fringe pattern, varying 
randomly between snapshots. This experimental result 
begs the old question [3| of how on earth two independent 
systems can exhibit knowledge of the relative phase of 
their wave functions. In the present context, this funda- 
mental problem was analysed in Ref. 3] and Ref. 4] : Us- 
ing various models and approaches, indeed it was shown 
that the behaviour of condensates with large number of 
particles is indistinguishable from the broken symmetry 
scenario where the condensate is described by a wave 
function with a definite, although unknown, phase. Since 
then, the problem has been attracting attention from 
many researchers (see 0, 0, 0, Q ) ^^'^ references therein) . 
The purpose of this paper is two-fold: First, we present 
a method, which we believe is most suitable for study- 
ing density fluctuations in macroscopic quantum systems, 
and second, we reanalyse the problem of interfering BEC 
clouds using the method. 

Interfering BEC clouds provide probably a unique ex- 



ample of a system with a large number of particles on 
which a measurement gives macroscopically distinguish- 
able results that are unpredictable. Indeed, the interfer- 
ence images registered in identical experiments differ in 
the fringe position, and the macroscopic phase extracted 
from the image is random from run to run. For this 
reason, the standard description of macroscopic systems, 
where fluctuations are small and expectation values are 
representative, is not sufficient, and the information we 
seek is not an expectation value, but the probabilities for 
realisations of possible density profiles of BEC clouds. 
Such detailed information is normally not of concern for 
macroscopic systems, but is precisely the information 
needed for discussing the outcome of a density measure- 
ment on interfering BEC's, viz. to extract a phase from 
an image, a global feature and an information which can 
not be extracted from local expectation values. 

Our goal is to obtain directly the statistics of instanta- 
neous measurements of particle distributions, the kind of 
information needed to describe snapshot experiments on 
many-body systems such as BEC's [l^. Given a desired 
spatial resolution lx,y,z in each of the three dimensions, 
we imagine a mesh in real space where the elementary 
volume is a Za; X Zj, X Iz box of volume A = Ixlyh- At the 
finest level, either a particle or none is detected in each 
such pixel. The result of such a position measurement 
of the particle locations in a many-body system corre- 
sponds thus to zeros and ones attributed to each such 
pixel. Coarse-graining, i.e., adding the observed parti- 
cle content of nearby cells, the extend of which specifies 
a bin, will attribute a corresponding integer number of 
particles to each bin. The result of a snapshot measure- 
ment of the positions of the particles is thus specified 
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by a set of integers^ {n\ = ni, n2, . . . , n^, . . . the num- 
ber of particles detected in the bins located at positions 
Xi,X2, . . .Xk, . . . respectively. The size of the bin (e.g. 
the pixel of a CCD camera) determines the resolution of 
the image of the particle distribution. 

Our aim is to find an efficient method for obtaining the 
probability distribution, P({n}), for the various strings 
{n} = ni,n2, ... of integer outcomes of the particle dis- 
tributions for a general quantum system consisting of N 
particles. 

To this end we employ the phase transformation of the 
wave function . . . , xjv) — > IJ\ ?A(xi, . . . , xat) where, 

similar to a gauge transformation, 

tlx ^(xi, . . . , x^) = e*(^(-i)+-+^(-"))V'(xi, . . . , xjv). 

(1) 

(A related phase tagging was employed in an investiga- 
tion of counting statistics developed for studying electric 
current noise in mesoscopic devices @.) Tagging in this 
fashion the particles in the state j?/'), we construct the 
following functional of the phase field 



$[A] = (V-IJ/aIV')- 



(2) 



The Fourier transform of <i>[A] with respect to the 
"gauge field" A(x) provides the probability distribution 
for strings, specifying the possible particle configurations. 

Indeed, let the real function A(x) be piecewise con- 
stant in each of the bins where the number of particles 
rik is counted, i.e., A(x) — Afc in the fc-th bin. The func- 
tional <i>[A] becomes a function, <i>({A}) of the variables 
Ai, A2, . . . = {A} The probability distribution, P{{n\), is 
then obtained as the multidimensional Fourier transform 
with respect to the A in each of the bins. 



P{{n})^ y" DA e-*^'="''^'=$({A}), 


where DA denotes the product of the differentials, 
dAi dA2 dAfe 



DA = 



2-K 2-K 



27r 



(3) 



(4) 



and the 27r-integration is performed with respect to each 
of the A's. 

To validate this construction, we note that in accor- 
dance with general rules of quantum mechanics the prob- 
ability P({n}), i.e., the joint probability that exactly 
particles are found in the fc-th bin, fc = 1, 2, etc, is given 
by the integral. 



P{{n})^ J dxi . . .dxA,|-0(xi, 



(5) 



where the region of integration, specified by the string 
{n} = 711,712,... is restricted by the set of conditions 
that for any k — 1,2, etc exactly Uk arguments of the 
wave function belong to the fc-th bin. 



The validity of the counting formula Eq. ([3]), where 
^({A}) is given by Eq. ^) with the piecewise constant 
A(x), is immediately seen since the Fourier transform in 
Eq. ^ has the effect that the unrestricted spatial in- 
tegration in Eq. ([2]) is narrowed to the region {n} as in 
Eq. This concludes our proof that the functional <i>[A] 
is the generating function for the particle dstributions. 

For future needs, we present Eq. ([3]) in a concise form 
as a functional integral with respect to the field A(x) 
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P[n] = J pAe-'/'^^"(^)-^(^) $[A] 



(6) 



where VX is understood as in Eq. and n(x) = 

Sfc ''^k ^(x — Xfc) specifies an arbitrary assignment of the 
bins filling or it may be viewed as a continuous density 
distribution ri(x). 

We have thus achieved the goal of expressing the statis- 
tics of particle locations, not as customarily directly in 
terms of the absolute square of the wave function, but 
in such a fashion that the definite outcome of particle 
content detected in each bin, Uk, enters explicitly as pa- 
rameters. The spatial resolution of the profile 77(x) is 
controlled by the smoothness of the field A(x) one chooses 
when evaluating the functional A-integration in Eq. 
The formula, Eq. ([3]), thus provides the complete count- 
ing statistics of single shot particle distributions for an 
arbitrary many-body system. As an application of the 
method, we consider the counting statistics of interfer- 
ing Bose-Einstein condensates. In this example, we re- 
strict ourselves to the limit of non-interacting particles, 
neglecting correlation effects related to the depletion of 
the condensates. 

First, we consider an isolated BEC cloud consisting of 
N bosons. For our purposes it will be convenient to build 
the Fock state of the cloud in the coherent state represen- 
tation. Let 0(x) denote the state into which condensation 
occurs, and consider the coherent state 



_ ^e«'/dx0(x)V;t(x)| 



|0), /dx|0(x)|2=iV (7) 
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obtained by operating with the boson field V'^(x) on the 
vacuum state |0).^1^| The Fock state \N) where N bosons 
occupy the single particle mode 0(x) is then obtained by 
projecting the coherent state, Eq. ([7]), 



(8) 



where cn = y^Nl/N^ is the normalisation constant. 

We then consider two initially separated clouds of 
bosons with particle numbers Ni and A^2, and where the 
particles occupy single particle states (/'i(x) and (f>2{x), 
respectively. The many-body wave function of the two 
spatially non-overlapping clouds is simply the product of 
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the wave functions of the isolated clouds, and the corre- 
sponding quantum state, |^) = |-/Vi) 11^^2)2, is thus ob- 
tained expressing the Fock state in terms of the coherent 
state, according to 



27r 



l^-) = C12 



where 



■09162 (x) 



j(eiAri+e2A'2)g/dx^eie2 WV-^WjO) 



'0l(x) 



^02(x) 



(9) 



(10) 



and C12 is the normalisation constant, C12 — cniCn^. 

One may get the wrong impression that in the quantum 
state VP, Eq. the bosons are in the superposition state 
^eie2(x) of Eq. pO)) . However, unlike genuine superpo- 
sitions, where a phase transformation like ~^ 
non-trivially changes the state, the same transformation 
of the many-body state in Eq. ^ amounts to an over- 
all phase factor multiplication which does not change the 
value of any observable. In fact, the atoms in the two 
non-overlapping clouds are not yet even aware that they 
are identical species. 

In the interference experiment, time evolution provides 
by the time of measurement, ioj overlap of the initially 
separated clouds due to their expansion. If the atom- 
atom interaction is neglected, evolution of the many- 
body state l^'(i)) amounts to the unitary evolution of 
the single particle states. Consequently, the expression 
Eq. ^ is valid for the state Y^it)) provided 0i^2(x) are 
substituted by the evolved single-particle wave functions 
(/)i.2(x, t). The initial states </>! and (j)2 are spatially sepa- 
rated, i.e., locally the states satisfy (/)i(x)02(x) — 0, and 
are therefore trivially orthogonal. It follows from the uni- 
tarity of evolution that the spatially overlapping states 
(/>i(x, t) and (/)2(x, i) remain orthogonal at any time. 

To simplify notation, we use Eq. ^ for the state 
1 4' (to)) of the clouds at the moment of measurement, t^, 
with the understanding that in all formulae below (/ii^2(x) 
stand actually for (/)i^2(x, to); the functions ipi and 02 are 
orthogonal and normalised to the particle numbers A'^i 
and N2 of the respective clouds. 

Having specified the state of the system, we are in the 
position to obtain the probability for any possible particle 
configuration [n] from Eq. ([6]) . To evaluate the functional 
$[A] = (5'|C/a|^') with I*) given by Eq. (O, we note that 
in terms of the field operator, the tagging of the particles 
by the phase field via the operation J7a|\E') is enforced by 
its transformation according to 



^(x) 



V'(x). 



Using the Baker-Hausdorff formula [10l| to calculate the 
matrix element (vI'|[/a|5') we arrive after straightforward 
calculations at the expression for the particle probability 



distribution 

X exp [/ dx (e*^Wv;'(x)V'e(x) -m(x)A(x))] 
where for brevity continuum notation is used, and 



06 (x) 



01 (X) 



02 (X) 



(11) 



(12) 



with Via Ni^/N, N =^ Ni + N2.[15] Eq. ^ refers 
only to configurations n(x) which are physical, i.e., 
/dxn(x) = N, and otherwise the probability is identi- 
cally zero as observed in the process of obtaining Eq. (fTT|) 
from Eq. ([6]). For brevity, in Eq. (jlip and below we omit 
the Kronecker symbol expressing this property. 

The formula Eq. pT|) . allows one to find the particle 
probability distribution for two Fock state clouds with 
arbitrary particle numbers Ni and N2. 

The functional A-integral is understood as the multi- 
dimensional one in Eq. ((3]). In our approximation, where 
the interaction is ignored, the integral factorises and can 
be easily evaluated. Appropriately choosing the size of 
the bins, any spatial resolution for the density profile 
n(x) can be considered. 

For illustration, we first look at the simplest case where 
each "cloud" consists of a single boson, A^i = 7V2 = 1, in 
the state 0i(x) and 02 (x), respectively. In this case, only 
two types of particle configurations exist: (i) both parti- 
cles are located in the bin around point x^; {n(x)}aa = 
Oi, O2, . . . , 2a, . . . where all particle numbers Uk are zero 
except for bin a where tIq = 2. (ii) One of the parti- 
cles is located in the bin near and the other one is 
in a different bin x;,; {n(x)}ah = Oi, O2, . . . , la, ■ • ■ lb • ■ •, 
i.e., rife = excepting Ua = rib = 1- Denote by Paa 
and Pab the probabilities for the corresponding particle 
configurations, and in order to single out the dependence 
on the bin volume present them as Paa = "Paa^^ and 
Pab = Vab^^ , assuming the bins small so that 0i,2(x) 
do not vary within a bin. To obtain the probabilities 
for these configurations, one integrates in Eq. pT|) with 
respect to the field A(x) assuming that the latter is piece- 
wise constant in each of the bins. After further integra- 
tion with respect to the phases 9 and 9' , one obtains 

Vaa = 2|0i(Xa)02(Xa)|', (13) 
Tab = |0l(Xa)02(xb) +0l(Xb)02(Xa)|^. (14) 

Here, the coefficient of 2 in the expression for Vaa repre- 
sents the well-known boson bunching effect [lOj]. The 
cross term, 23fi0i(xa)0t (xb)02(xb)02(xa), in Tab de- 
scribes the expected Hanbury Brown and Twiss effect 
[Toj . Also, one observes that 'Pab\b^a — 2Vaa', this is a 
necessary property to guarantee the correct dependence 
of Paa on the choice of the mesh size. 

Next we turn to the case of interest, two interfering 
macroscopic Bose-Einstein condensates where the num- 
ber of particles in each cloud is large, A^i^2 ^ 1- The main 
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contribution to the phase integral in Eq. (fTTj) comes in 
that case from the region 6' m 6, and the particle proba- 
bility distribution becomes 



Z7T 



(15) 



where 



Vg[n] = V2^JVX e/dx((e'^W-l)p«(x)-»n(x)A(x)) 

with 

pe(x) = |e*Vi(x) + ^2(x)|' . 



(16) 
(17) 



Within the approximation of non-interacting bosons 
the formula Eq. (|15p provides the complete counting 
statistics of snapshot particle distribution measurements 
of two macroscopic interfering Bose-Einstein condensates 
in Fock states, i.e., a state where the particle numbers 
iVi,2 of the condensates are large and separately fixed 
and their phases are completely indefinite. 

To get insight into the content of formula Eq. (fT5|) , we 
note that the distribution P[n\ is built by a family of 
partial (normalised) distributions Vein]. Obviously, the 
part of the configuration [n] -space where the probability 
P[n\ is appreciable, is spanned by the corresponding sub- 
spaces generated by the Ve[nys, < 9 < 2tt, and we first 
discuss the partial contributions Ve [n] . One can check by 
comparison that Ve [n] coincides with the distribution of 
the outcomes of N independent identical double-slit-type 
experiments where the position x of a single particle in 
the superposition state V'e(x) = e'^0i(x) + (j)2{x) is mea- 
sured. At the individual pixel level, the distribution of 
the number of particles n^, registered in the fc'th bin is 
Poissonian reflecting the intrinsic randomness at the level 
of quantum mechanics. When typical counts n^'s are 
large, the shot noise fluctuations in the n^-'s are relatively 
small. Asymptotically, at N oo, Ve[n{x.)] thus selects 
the configurations where the density n(x) equals (up to 
shot noise) the density expectation value with respect 
to the wave function building /5e(x) in Eq. (jl7p . Conse- 
quently, the distribution P[n] becomes concentrated on 
the (closed) line corresponding to [n] — pe with 9 running 
from to 2tt. Excepting exponentially rare events, the 
outcome of any experiment will correspond to a point in 
the close vicinity of this line, and consequently, an exper- 
iment will always show the interference picture as if the 
particles are condensed in a state with a certain definite 
phase difference 0, |3| • 

It is well understood [1, 0, d, d, 0| that each experi- 
mental run results in a fringe image which corresponds to 
a well defined phase difference 9, but its value is random 
when the experiment is repeated under identical condi- 
tions. Having at hand the full counting statistics of the 
interference image, we are able to go further into such an 
analysis and discuss the probability, p{9), of the occur- 
rence of the images corresponding to the phase 9. Unlike 



previous authors, we do not consider it a priori obvi- 
ous that any value of 9 is equally probable as we are 
not aware of any physical symmetry that demands states 
with different phase difference 9 to be equivalent. On the 
contrary, we observe that the entropy of the partial distri- 
bution Ve is ^-dependent. i.e., the statistical weight 
of configurations to which a certain value of 9 can be 
assigned varies with 9. Besides, there is no unique and 
unequivocal algorithm of the assignment, and the ques- 
tion warrants to be addressed. In general terms, this 
is a standard problem of mathematical statistics where 
one determines the parameters of a distribution from a 
sample data [ll[. 

Specified a procedure Q[n] that assigns the phase dif- 
ference 0[„] = Q[n\ to a configuration [n and given 
the occurrence frequency of the configurations P[n\ in 
Eq. (fTSl) . the probability distribution p(0) is found as 



where 



PiO) = / ^We-i^) 



(18) 



(19) 



is a function of 9, normalised to unity, and with 9' as pa- 
rameter. For an unbiased estimator Q[n\ 11] , for which 



J2Q[n]Ve[n] ^ 9 , 



(20) 



the function W is strongly peaked at |6' — 6*'! ^ ae', where 
the width of the peak, ae', can be estimated from the 
Cramer- Rao bound ^11] , 



dxp0(x) 



din pg(x) 
09 



(21) 



The physical meaning of the variance ag is that of the 
confidence ( "the error bars" ) with which one is able to es- 
timate the phase difference from a noisy image. By virtue 
of the normalization, pe is proportional to the number of 
particles in the clouds so that the variance ae becomes 
small in the macroscopic limit. Note that the variance 
varies with 9 in agreement with the above discussion of 
the dependence of the statistical weight on the phase. 
In the large- limit, the function We' (9) is a Gaussian, 



We' (9) 



1 



V^7 



-[e-6'fl1al, 



the normalisation property with respect to 9 follows from 
Eq. (fT9|) . Although the variance agi is a function of the 
phase difference 9' , in the integral with respect to 9' , 
Eq. (|18p. one can safely replace ae' — > ae, after which 
the integral becomes a constant, and p(6') — l/2n is in- 
dependent of 9. Thus, we see that the statistical inequiv- 
alence of the partial distributions Ve is compatible with 
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the statement that an image of two overlapping conden- 
sates can be characterised by a phase difference 9, the 
value of which in repeated experiments under identical 
conditions is equally probable in the interval < < 27r. 

Our approach allows us to evaluate the precision 59 
with which the phase difference may be meaningfully as- 
signed to a noisy fringe image: Using Eq. P?|) . the width 
of the "error bars" given by <Tg in Eq. (I2ip . can be ex- 
pressed via the wave functions (/)i.2([x) of the interfering 
clouds: 

-2_1 l0iPl02p-SR((e'Vi0^)^) 

' "27 ^''|0iP + |0,|2 + 25R(eiVi02)' ^ ^ 

For a simple illustration, we consider clouds which are 
initially in the ground states of ID harmonic traps cen- 
tered at x = ±A. After a free expansion during 
time r, the wave functions of the clouds are 01,2 (a;) = 

y/Ni^2/{aT\/^) exp [—{x ± A)^/ {2aar)] where Ur ~ a + 
ihT/{ma), m and a being the particle mass and initial 
width of the clouds, respectively. If the clouds overlap is 
small, i.e., \ar\ < A^/a the uncertainty in the phase 69 
is estimated as 69 ^ ae ^ exp [A"^ / {2a\ar\)] /{NiN2y/^. 
As no surprise, the uncertainty of the phase parameter 
69 oc {NiN2)^^^^ is small provided A^i,2 is large enough. 

In this paper, we have presented a method which allows 
one to evaluate the probability of a local density config- 
uration, [n] — n(x), for a general many-body quantum 
state. The configurations are specified by a string of in- 
tegers {rifc}, fc = 1, 2, etc , where the number at the fc-th 
position is the number of particles located in the cell cen- 
tred at position x^. The cells cover the relevant part of 
space, and the size of the cells defines the resolution of the 
density distribution. The probability of a density profile 
is expressed via the functional integral over the phase 
field A(x), Eq. © or Eq. ©, and the functional $[A(x)] 
defined by Eqs. H]) and As usual in the probabil- 
ity generating function approach, the n-point correlation 
function can be obtained by taking the functional deriva- 
tives of the functional $[A(x)] with respect to A(x) 
at the corresponding n values of the co-ordinate x. We 
have considered the case of only one species of particles, 
the generalisation to many species, bosons or fermions, 
is straightforward. 



As an example, we have applied the method to the de- 
scription of the interference experiment with ultra-cold 
atoms, where the quantum state is formed by two over- 
lapping but initially separated clouds of Bose-Einstein 
condensates. We describe the condensates in the approx- 
imation of non-interacting bosons and, in particular, we 
ignore the condensate depletion. At the level of self- 
consistent description, the Gross-Pitaevskii theory, the 
interaction only modifies the evolution of the interfer- 
ing waves (/)i,2(x) as discussed in |T2|. Beyond the self- 
consistent approach, the most important effect of inter- 
action is nontrivial spatial correlations in the noise of in- 
terfering BEG images (see Ref. [8i] and references therein) . 

The density profile of the interfering matter waves gen- 
erated by the two independent sources shows strong fiuc- 
tuations, and the standard description of macroscopic 
systems in terms of the expectation value of the den- 
sity and its low order, e.g. 2-point, correlators becomes 
insufficient. The "hidden" variable of interest is the pa- 
rameter 9 which corresponds to the phase difference of 
the matter waves in the broken symmetry BEC-state, 
and this parameter can be extracted by analysing the 
density profile, i.e., the CCD image registered in a snap- 
shot experiment rather than from the expectation value 
found by averaging many such measurements. We have 
evaluated the probability for arbitrary density profile re- 
alisations as well as the probability distribution p{9) for 
the phase difference corresponding to these profiles. Sup- 
porting expectations of previous authors d, 0, d, d, 0] , 
we have found that the phase difference 9 is distributed 
uniformly for an unbiased procedure of assignment of 9 
to the observed images. 

In conclusion, we believe that the approach developed 
in this paper offers an efficient tool for a detailed analysis 
of the problem of interfering matter waves and gives new 
insight into this fundamental physical problem. 
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